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ABSTRACT 

^ ■ 
O ; 

^ ! We analyze the statistical properties and dynamical implications of galaxy dis- 

5_i ! tributions in phase space for samples selected from the 2MASS Extended Source 

1^ ! Catalog. The galaxy distribution is decomposed into modes 5(k, x) which de- 

^ \ scribe the number density perturbations of galaxies in phase space cell given by 

^ ! scale band k to k + Ak and spatial range x to x + Ax. In the nonlinear regime, 

! 5(k, x) is highly non-Gaussian. We find, however, that the correlations between 

^ ! 5(k, x) and 5(k', x') are always very weak if the spatial ranges (x, x+ Ax) and (x', 

OO \ x' + Ax') don't overlap. This feature is due to the fact that the spatial locality of 

5^ \ the initial perturbations is memorized during hierarchical clustering. The highly 

2 \ spatial locality of the 2MASS galaxy correlations is a strong evidence for the 

Tj- \ initial perturbations of the cosmic mass field being spatially localized, and there- 

^ \ fore, consistent with a Gaussian initial perturbations on scales as small as about 

'"q^I 0.1 h~^ Mpc. Moreover, the 2MASS galaxy spatial locality indicates that the 

^ \ relationship between density perturbations of galaxies and the underlying dark 

^ \ matter should be localized in phase space. That is, for a structure consisting 

^ \ of perturbations on scales from /c to /c + A/c, the nonlocal range in the relation 

^ \ between galaxies and dark matter should not be larger than |Ax| = 2'7r/|Ak|. 

\ The stochasticity and nonlocality of the bias relation between galaxies and dark 

^ \ matter fields should be no more than the allowed range given by the uncertainty 

■ ■ ' relation |Ax||Ak| = 27r. 

Subject headings: cosmology: theory - large-scale structure of the universe 



1. Introduction 

The large scale structure of the universe developed from initial mass density and ve- 
locity fluctuations through gravitational instability. Much of the information of the initial 

^Center for Astrophysics, University of Science and Technology of China, Hefei, Anhui 230026, P.R.China 
^Department of Physics, University of Arizona, Tucson, AZ 85721 



- 2 - 



perturbations is "forgotten" during the gravitational nonlinear evolution. Yet, some features 
of the initial perturbations are imprinted in the cosmic mass and velocity fields at present. 
Measuring these memorized features is crucial in studying the initial states of the universe. 
A well known example is the self-similarity of gravitational clustering, which imprints the 
initial power spectrum index and is detectable with the scaling behavior of correlation func- 
tions of the present mass field (e.g. Peebles 1980). In this paper, we study the spatial locality 
of the perturbed mass field in phase space, which is also a feature to be memorized in mass 
field today. 

The physics of spatial locality of correlation function in phase space can be illustrated 
with a Gaussian mass field p(x) described by 

m)s{-k')) = p{k)6^,^„ (1) 

where 5(k) is the Fourier counterpart of the density contrast 5(x) = [p(x) — P{k) 
the power spectrum of the mass field, and the Kronecker delta function. Eq.(l) says 
that the Fourier modes with different wavevcctor k are uncorrelated, or the correlation is 
localized in A;-spacc. This is because the phase of the Fourier modes 6(k) is random. On the 
other hand, the correlation function of density perturbations in physical (x) space generally 
is non-local. The two-point correlation function (5(x)5(x')) has non-zero correlation length 
when the Fourier power spectrum P{k) of cq.(l) is fc- dependent. 

In a phase-space description, the mass field is decomposed into modes 5(k, x), the 
perturbations in the wavevcctor (scale) from k to k + Ak and physical range x to x + Ax. 
The volume of the phase space cell referring to a mode is given by the uncertainty relation 
I Ax 1 1 Ak I = 27r. The correlation function of a Gaussian field generally is localized regardless 
whether the Fourier power spectrum is colored. That is, 

(<5(k,x)5(k',x'))oc5j,,5£,,. (2) 

The reason for eq.(2) is straightforward. First, the perturbations (5(k, x) and 5(k',x') are, 
respectively, given by linear superposition of the Fourier modes in different wavebands (k, 
k -|- Ak) and (k', k' -|- Ak'). For a Gaussian field, the Fourier modes in different wave band 
are uncorrelated in general [eq.(l)]. This gives rise to the factor of 5^^^, of eq.(2). Second, the 
phases of the Fourier modes of Gaussian field are random. For a superposition of the Fourier 
modes in the band (k, k -|- Ak) with random phase, the coherent length of the phase of 
5(k, x) is not larger than |Ax| ~ 27r/|Ak|. On the other hand, the spatial distance between 
two cells in the (k — x) space is at least |Ax| ~ 27r/|Ak|. This yields the factor 5^^, of eq.(2). 
Therefore, the phase space correlation function eq.(2) for Gaussian field generally is local. 

Non-linear evolution via gravitational clustering will lead to a non-Gaussian field that 
will deviate from the locahty of eq.(2), even when the initial field is Gaussian. However, 
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turbulence studies have found that if 1.) the initial perturbations are spatially localized, 
and 2.) the random fields evolve via a self-similar hierarchical cascade process, the phase 
space correlation function of the evolved field will still be spatially localized (Grcincr, Lipa, 
& Carruthers, 1995, Greiner et al. 1996). The perturbations of modes 5(k, x) and 5(k',x') 
with X 7^ x' stay statistically uncorrelated or only very weakly correlated during the self- 
similar hierarchical cascade evolution. That is, the factor 6^^, of eq. (2) is memorized during 
the dynamical evolution. 

Phenomenological models that mimic the hierarchical clustering of the cosmic mass field 
have similar mathematical structures as the hierarchical cascade models of turbulence. For 
instance, the fractal hierarchy clustering model of Soneira & Peebles (1977) is the same as 
the P model of turbulence (e.g. Frisch, 1995). The block model of Cole & Kaiser (1988) is a 
special case of the multifractal cascade model (Meneveau & Sreenivasan 1987, Pando et al 
1998). Hence, these models should also memorize the spatial locahty. Recently, the spatial 
locality has been studied with more realistic dynamical models of gravitational clustering. 
First, if the weakly nonlinear mass field is given by the Zcldovich approximation, the field is 
found to be spatially localized if the initial perturbations are Gaussian (Pando, Feng & Fang 
2001). More recently, this result has been extended to fully nonlinear regime (Feng & Fang 
2004). Using the halo model of the large scale structure (e.g. Cooray & Sheth 2002, and 
references therein), it has been shown that the evolved mass field is approximately spatially 
localized if the initial perturbations are Gaussian. Although gravitational coupling is long- 
term, the spatial locality in the phase space is not disturbed by the non-linear evolution of 
cosmic mass field. This property essentially is due to the self-similarity of the hierarchical 
clustering. This result has been tested with high resolution N-body simulation samples (Feng 
& Fang 2004). 

In this paper, we investigate the spatial locality of phase space correlations with real 
sample - the galaxies selected from the 2MASS Extended Source Catalog (XSC, Jarrett 
et al. 2000). Our motivation is two-fold. First, the spatial locality provides a test of the 
Gaussianity of the initial density perturbations on small scales. Although many tests on the 
Gaussianity of the initial perturbations have been done with the temperature fluctuations 
of the Cosmic Microwave Background Radiation (e.g. Komatsu, et al. 2003, Pando, Valls- 
Gabaud & Fang, 1998), the comoving scales of these tests are not less than a few Mpc. The 
spatial locality of the 2MASS samples can test the Gaussianity to scale as small as about 
0.04 h~^ Mpc. Second, the distribution of galaxies is biased from the mass field of dark 
matter. Some bias models assume that the relation between the distribution of galaxies 
and underlying mass field is stochastic and nonlocal (e.g. Dekel & Lahav, 1999). This 
mechanism will lead to nonlocality of the galaxy correlation function in phase space, even 
when the underlying dark matter mass field is spatially localized. Therefore, the spatial 
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locality should be effective in testing the stochasticity and nonlocahty of bias models. 

The outline of this paper is as follows. §2 presents the statistics and dynamics of the 
spatial locahty in phase space. §3 describes the basic properties of the 2MASS galaxy samples 
with a space-scale decomposition. The results of the spatial locahty of the 2MASS galaxy 
distribution are presented in §4. Finally, the conclusions and discussions are be given in §5. 

2. Statistics and dynamics of spatial locality in phase space 

2.1. Variables of the mciss field in phase space 

The cosmic mass density contrast field 5 in x-space is given by 5(x) = (x|5), in k-space 
by (5(k) = (k|(5). These two descriptions are equivalent, and the set of the bases |x) or |k) 
are complete and orthogonal. S{x) and 5(k) correspond to, respectively, the coordinate- and 
momentum-representations of quantum mechanical wavef unctions. To study the statistical 
and dynamical behavior of the mass field in k — x phase space, we should describe the field 
5 in the Wigner representation. Therefore, the phase space analysis of 5 can be performed 
with the discrete wavelet transform (DWT). The Wigner wavefunction is a prototype of the 



For the DWT analysis, we use the same notation as Fang & Feng (2000) or Fang & 
Thews (1998). In the DWT scheme, there are two sets of base given by 1.) scaling functions 
0j^i(x) = (x|j, l)s and 2.) wavelets ipj,i{'x) — (x|j, where j — 0, 1, .. and I — 0, 1...2^ — 1. 
In 1-D space with size L, the scaling function (f)j,i{x) is localized in physical space < 
X < {I + l)L/2^ , while wavelet ipj^i is localized in phase space cell lL/2^ < x < (l + l)L/2^ 
and Ti2^ / L < \k\ < (3/2)27r2''/L. The DWT base in 3-D space is given by a direct product 
of the 1-D base, i.e. |j,l)s = /i)s|j2, /2)s|j3, ^3). and |j,l)^ = |ji, /i)^|j2, /2)«;|j3, ^3)«>- If 
a sample is in a cubic box of < Xj < L, i = 1, 2, 3 and volume V = L^, the bases 
are localized in the spatial range liL/2^^ < Xi < {U + l)L/2^^ and wavenumber range 
7r2^*/L < \ki\ < {3/2)2tt2^^/L. The simplest scahng and wavelet are given by 



DWT. 





w 





otherwise 



This is the so-called Haar wavelet. However, we will not use the Haar wavelet in the numerical 
calculation of §3 and 4, because the discontinuity of the Haar function made it not to be 
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very well localized in the Fourier space (/c-space). We will apply the wavelet of Daubechies 
4 (Daubechies, 1992), which has much better behavior in k space. The scaling and wavelet 
of the Daubechies 4 in x- and /c-spaces can be found in the Numerical Recipes (Press et al 
1992) and Yang et al (2001). The properties of the scaling function and wavelet discussed 
in this section are generally available for either the Haar, Daubechies or other DWTs. 

The scaling functions are orthogonormal with respect to index 1 as 

= y 0j,i(x)0j,i/(x)(ix = ^1^,. (3) 

Scahng function 0j,i(x) actually is a window function for the spatial range < Xi < 

{li + l)L/2^\ The scaling function coefficient (SFC) of a density field p is defined by \ = 
s(j, = / p(x)0j i(x)(ix, and therefore, the mean density in 3-D spatial range of < 
Xi < {k + 1)1/ 2^^ is 

where |1) is a uniform field with density equal to unity. The density field can be expressed 
by the SFCs as 

p(x)=p(J)(x) + 0(>j), (5) 

where k runs 0, 1,...2-'*— 1 and 0{> j) means all fiuctuations on scales less than (L^/2-'i+-'2+-?'3)i/^, 
and 

P^^Hx) = E^jAi(x)- (6) 
1 

Therefore, p^j^(x) is the density field smoothed on scale (L^/2-'i+-'2+J3^i/3_ ^j^g SFCs ej^i 
and are similar to the mass field variables given by count-in-cell. Thus, the two-mode 
correlation will display the similar features as ordinary two-point correlation func- 

tion, which generally is not localized for Gaussian fields with a /c-dependent Fourier power 
spectrum P{k). 

The wavelets '0j,i(x) are orthogonal with respect to both indexes j and 1 

= y"V'j,i(x)V'jM'(x)rfx = (7) 

The wavelet function coefficient (WFC) of the density contrast field 6 is defined by 

ej,i= = J (5(x)V'j,i(x)dx. (8) 

Since the set of the wavelet bases |j, is complete, the density contrast field can be expressed 
as 
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where each ji runs 0, 1, 2... and k runs 0, 1,...2^^ — 1. Therefore, the WFCs ej^i can be used as 

the variables of the mass field S. cj i is the fluctuations of the density field around scales k = 
27rn/L with n = {2^\2^\ 2^3) and at the physical area 1 with size Ax = {L/2^\L/2^^, L/2^-'). 
In each dimension i, iji is given by a superposition of fluctuations in the waveband ki±Aki/2, 
where ki = 2n2^^/L, and Aki = 2tt/ Axi = ki. Therefore, ej,i can play the role as the variables 
5(k, x) of eq.(2). We have then 

(^j,i^jM') = PiS^yS^v (10) 

As we have emphasized, the locality of the factors S^y and S(\, in eq.(lO) generally are very 
good approximation for Gaussian fleld. The factor Pj in eq.(lO) is the DWT power spectrum 
of the mass density perturbations (Fang & Feng 2000). Pj is not dependent on 1, as cosmic 
mass fleld is assumed to be randomly uniform. 

The spatial locality 5/^, arises from the phase decoherence given by the superposition 
of the Fourier modes in the band k to k + Ak. This point can easily be seen from the 
orthogonal relation eq.(7) by taking j = j'. For the 1-D case we have 

^5' = J ipjAx)tp3,i'{x)dx (11) 

^ oo 

n=—oo 
^ oo 



2J 



. J2 V^'(n/2^>-^2'^'^('-'')/2', 



n= 



where k — 27m/ {L/2^), and 'ipj,i{k) is the Fourier counterpart of ■0j,«(^)- The last step of 
deriving eq.(ll) used the following relation (Fang & Thews 1998) 

iy(^)=(^)'^'v^(r^/2^)e-'^'^'^'/^', (12) 

where i){z) is the Fourier transform of the basic wavelet. It is non-zero only within z — l±l/2 
and -1 ± 1/2. Therefore ^(n/2-') is non-zero only for 2-'(l - 1/2) < |n| < 2^{l + 1/2). Since 
the factor %lP'{n/2^) is always positive, the spatial locality factor 5/^, of eq.(ll) basically 
comes from the average of the phase factor e~*^'^"'('~'')/^^ over n. For (/ — /') ^ 0, the phase 
27rn(/ — l')/2^ is uniform distributed in range (l/2)27r — (3/2)27r. Therefore, the phase is 
similar to a randomly uniform distributed. Thus, for a function /(x), we have approximately 



/ 



^ oo 

^,,{x)f{x)^,,r{x)dx ^-Y. fm\nl2^)e-'^-<'~'^'^' cx (13) 



2J 

n=— oo 



if f{k), the Fourier counterpart of f{x), does not change the average of the phase factor 

^-i2^n{l-l')/2i ^jjg j.^^gg 2^(1 _ 1/2) < \n\ < 2^(1 + 1/2). 
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2.2. Statistical criterion's for spatial locality 

The mass field evolution under a self-similarly hierarchical process will not violate the 
initial spatial locality (Feng & Fang 2004). The correlations between DWT modes at different 
physical positions are always substantially less than that at the same position. In other 
words, gravitational evolution causes a strong coupling between the phase space modes (j, 1) 
and (j', 1') with different scales j and j', but with spatial range of 1 and 1' overlapped. The 
coupling between modes at different physical range is very weak. The factor in the initial 
Gaussian field [eq.(lO)] is memorized during the nonlinear evolution. 

The spatial locality can be measured by 

mA =/i 1=1' (14) 

where (...) is average over the ensemble of variables ij i with /j = 0...2'* — 1 for a given j. 
Considering the mass field p is randomly uniform, we have (|ej,ip) = Therefore, 
eq.(14) can also be rewritten as 

|/c;x/c?xVj,i(x)e(x-xOV;j,i,(xO| 

|/cix/cix'^j,(x)e(x-x')^j,(x')| ' ^ ■ ^ ^ 

where ^(x — x') is the two-point correlation function of mass field p. Eq.(15) is a criterion 
of the spatial locality with second order statistics. 

One can construct the spatial locality with higher order correlations of the ej^i. For a 
ip + q) order statistical criterion, we use statistics as 

where p and q can be any even number. Obviously, for Gaussian fields, (ejie? j,) = (ej'i)(ej j,) 
at 1 ^ r, and therefore, C['« = 1 at 1 ^ 1'. If Cf^ ^ 1 at 1 ^ 1', the field is nonlocally 
correlated. Therefore, the spatial local evolution is given by 

Cf'-l ff 1^1'. (17) 

Eqs.(14) and (17) are the basic statistical criterion's used for testing the spatial locality. 



2.3. Bias model of galcLxies and spatial locality of correlation 



The number density of galaxies Pg(x) are not simply proportional to the mass density 
p(x) of dark matter. It is biased. The simplest linear model of galaxies bias is given 
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by 5g{x) — 65 (x), where Sg{x.) is the number density contrast of galaxies, and b the bias 
parameter. This model implies that the relation between 5g(x) and 5(x) is deterministic and 
localized. It has been argued that the locality assumption is not trivial, and consequently, 
stochastic and nonlinear bias models have been proposed. In these models, the galaxy field 
^^(x) at X is not determined only by dark matter field 5(x) at the same point x. The 
relation between 5p(x) and 5(x) is stochastic and non-local. In these models, the two-point 
correlation function of galaxies ^^(x — x') is related to the two-point correlation function 
^(x — x') of dark matter by (Matsubara 1999) 

^g(x — x') = J K{y)K{y')^ (x — x' -|- y — y')dydy' -\- higher order terms, (18) 

where K{y) describes the nonlocal nature of the bias. A galaxy at x is dependent on the 
dark matter field 5(x -|- y) with probability proportional to K{y). If K{y) is non-zero in the 
range |y| < D, the relation between galaxies and dark matter field is nonlocal with size D. 
Thus, the galaxy correlation of phase space modes becomes 

(ej,iej,r) = J ^j,i{x)^g{x - x')'4)j^i>{x')dxdx' (19) 

oo 

J2 K{k)P{k)ij\n,/2^')i>\n2/2^') 



ni,n2,n3=— oo 

where i^(k) is the Fourier counterpart of i^(y). -^(k) is non-zero only in the band |k| < 
27r/L'. 

Thus, if D is large than L/2'', the non-zero range of A'(k) is limited in < L2^^ / D 
or ni/2^^ < 1. In this narrow range of n^, the average of the phase factor e^*^™i(^'~'i)/^^* 
cannot be zero when |/j — /-I > 1. Therefore, even when the underlying field of dark matter 
is spatially localized as eq.(15), the galaxy clustering with nonlocal bias will not be spatially 
localized. Hence, if the spatial locahty holds for a given scale j, the nonlocal size of galaxy 
bias should not be larger than L/2^ , i.e. the nonlocal scale D will not be larger than that 
given by the uncertainty relation Ax — 2T:/k. 



2.4. 2-D distributions 

Our sample of the 2MASS galaxies is 2-D. With the DWT, the projection of 3-D field 
5(xi,X2,X3) to 2-D 5{xi,X2) is given by 

5{xi,X2)^ s{j3,0\S) ^ J S{Xl,X2,X3)(f)j^^o{Xs)dXs, (20) 
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where Xi and X2 are coordinates on sky plane, while X3 stands for the spatial coordinate in 
redshift direction. In eq.(20), 0^3^0(2^3) is a scaling function on scale Hence, if 

is equal to the depth of the sample in the redshift direction, eq.(20) is a projection of a 3-D 
sample 6{xi,X2,X3) with depth onto a 2-D distribution. The DWT variables in 2-D 

analysis is then 

e?f = J S{xi,X2)tlJji,ii{xi)tjjj2^i2{x2)dxidx2 (21) 

= J S{xi,X2,X3)'lpj^,i^{xi)'ljjj^,i^{x2)(l)jsfi{x3),dXidX2dX3. 

Eq.(21) shows e|f is a sampling of 5(x) by basis \ji, li)w\j2, ^2)«;|j3, 0)s, which is mixed from 
wavelets \ji,h)w, \32-,h)w and scaling |j3,0)s. 

Since the \ji,li)w and \j2,l2}w are orthogonal with |j3,0)s, the DWT is effective in 
distinguishing the redshift direction. One can repeat the analysis of Feng & Fang (2004) 
to show that the correlation of the variables ejf will be spatially localized with respect to 
the position indexes li or I2 regardless the projection of the redshift direction. The reason 
is simple. The scaling function base |j3,0)s counts all galaxies in the redshift direction. 
However, the wavelet bases and |j2,^2)s measure the difference heiween the densities 

of two neighboring cells on scale ji + 1 and ^2 + 1, and therefore, the background of the two 
neighbor cells are cancelled. The background may cause uncertainty of the shot noise, but 
it does not contribute to the difference e|f . Thus, if a 3-D mass field is spatially locahzed, 
the 2-D variables ^^ should satisfy eqs.(14) and (17). We do not write the superscript 2D 
below to simplify the notation. 



3. Samples for the spatial locality analysis 

3.1. Galaxies selected from 2MASS-XSC 

The galaxy samples used for the spatial locality analysis are selected from the 2MASS 
extended source catalog (XSC, Jarrett et al. 2000), which covers almost the entire sky 
at wavelength between 1 and 2 yum. To select galaxies, we use the indicator K_m_k20fe, 
which measures the magnitude inside a elliptical isophote with surface brightness of 20 mag 
arcsec"^ in i^^^-band (from then on, we infer K_m_k20fe as Kg). There are approximately 
1.6 miUion extended objects with Kg < 14.3. Most of the XSC sources at \b\ > 20° are 
galaxies (> 98%). The contamination mainly is from stars. The rehabihty of separating 
stars from extended sources is 95% at |6| > 10°, but drops rapidly to < 65% at |6| > 5°. 
To avoid this contaminant, we use a latitude cut of |6| > 10°. We also removed a small 
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number of bright {Ks < 9) sources by the parameters of the XSC confusion flag (ccJiag) 
and visual verification score for source (vc). They are identified as non-extended sources 
including artifacts. Moreover, to eliminate duplicate sources and have a uniform sample, we 
use the parameters use_src = 1 and dup_src = ^. 

To select the range of Kg, we use the standard log — log 5" test to examine the com- 
pleteness of the sample. The number counts can be approximated by a power-law as 

V oc 10'^"^. (22) 
am 

The XSC sources with |6| > 30° and 12 < Kg < 13.7 are believed to be galaxies with 99% 
rehability (Mailer et al. 2003). For this sample, the index k is found to be 0.641 ± 0.006. If 
considering this k to be the standard, the completeness of a sample can be estimated by the 
deviation of {dN/ dm) sample from the standard (Afshordi, Loh & Strauss 2003), i.e. 

{dN / dm) standard 

C{m) = , , , . (23) 

{dN/ dm) sample 

Figure 1 shows the number counts dN/dm and completeness C{m) for sample (1) (the 
standard) of galaxies with \b\ > 30° and 12 < Kg < 13.7, and sample (2) with \b\ > 10° 
region when 9 < Ks < 14. The figure shows that in the range 11 < Ks < 13.7, sample (2) 
has the same number counts and completeness as sample (1). The completeness C{m) of the 
sample (2) is obviously larger than sample (1) when Ks < 10.0. This indicates the catalog 
to be contaminated towards the bright end. C{m) drops below 0.9 when Kg > 14.0. Thus 
we use a cut of 10.0 < Ks < 14.0 to ensure our sample to be complete greater than 90%. 
This sample contains 987,125 galaxies. 

To carry out a 2-D DWT analysis, we first to take an equal-area projection with the 
Lambert azimuthal algorithm: 



xi = i?V2-2|sin6|cos;, (24) 
X2 — R\l2 — 2\ sin h\ sin 

where R is a relative scale factor, h is the Galactic latitude and / is the Galactic longitude. 
This hemisphere scheme projects the whole sky into two circular plane, northern and southern 
sky. From each circular plane, we select 14 squares, each of which has an angular size of 
about 28.28° x 28.28° (800 square degrees). The angular size labelled by j is 28.28/2^'. 



^The notations of the 2MASS parameters used in this paragraph are from the hst shown in the 2MASS 
Web site http://www.ipac.caltech.edu/2mass/releases/allsky/doc. 
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These squares do not overlap each other to guarantee the independence of statistics. Each 
square contains ~ 28,000 to 41,000 galaxies, and 35,000 in average. The 28 (28.28° x 28.28°) 
squares are our samples for statistics. Since XSC galaxies are resolved to lO", our analysis 
can reach to angular scale of 0.01 degree. For each square, we also produce random samples 
by randomizing the coordinate {xi,X2) of galaxies. 



The 2-D number density distribution of N galaxies with coordinate (x", X2) {n — 1...N) 
is given by 



where 6 is the Dirac delta function. We first calculate the correlation between modes of 
the density variable eq.(4), The variable pj^i is the mean number density in a 2- 

D spatial range < Xi < + l)L/2^\ i — 1,2. As have been discussed in §2.1, the 

correlation function is similar to the ordinary two-point angular correlation function 

w{9), where 9 is angular distance. A difference between w{9) and (pj,iPj,i') is that w{9) is 
calculated by counting the pair within a given angular distance, while the DWT algorithm is 
based on variable pj i corresponding to cell (j, 1) in phase space. It has been pointed out that 
the method of counting pairs actually is pair- weighted (Strauss, Ostriker & Ccn 1998). That 
is, a cell (j, 1) at the dense regions will be counted more than one time in the pair- weighted 
statistics, while empty cells contribute nothing. On the other hand, in the DWT analysis, 
each phase space cell (j,l) supports only one variable pj^i regardless the cell is dense, or 
empty of galaxies. 

Figure 2 presents the correlation function for mode j = (j, j) and j = 8, 9, 10 

and 11, corresponding, respectively, to smooth the sample on angular scale 28.3/2-' in unit 
of angular degrees. The angular distance in Fig. 2 is given by ^ = 28.3|1 — l'\/2^. The error 
bar is one sigma among the samples of 28 correlation functions given by the 28 data squares. 
As expected that the ^-dependencies of the correlation function show the typical power law 
in the angular range from ~ 0.01 to 1 degree. For j — 8, 9, the power law breaks at ~ 1 
degree, while for j = 11 correlation function is of a power law till scale ~ 2 degree. The best 
fitting results for a power law = A9~^ are shown by the solid lines in Fig. 2. The 

amplitude A and index /3 are listed in Table 1. 

Table 1 shows a systematic decrease of the index /? with the increase of j. Since the 
scaling functions are low pass filters, the SFCs at j smooth out all fluctuations on angular 
scales less than 28.3/2-^. Therefore, the j-dependence of /3 indicates that the clustering of 
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Two-Point correlation functions 




(25) 



n 
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Table 1: Power law fitting of correlation function 



j 


amplitude A 


index (3 


8 


0.154 ± 0.063 


1.020 ± 0.033 


9 


0.183 ± 0.002 


0.845 ± 0.017 


10 


0.190 ± 0.002 


0.819 ± 0.012 


11 


0.200 ± 0.004 


0.800 ± 0.008 



the 2MASS galaxies is stronger (and f3 is smaller) at smaller scales (or larger j). This 
result is consistent with the general picture that smaller scale clustering enters the nonlinear 
regime earlier. The density field of the 2MASS galaxies on small scales has significantly 
undergone nonlinear evolution. The values of P in Table 1 are httle larger than that given 
by the ordinary two-point angular correlation function (e.g. Mailer et al 2003). This is due 
probably to the fact that the pair- weighted algorithm gives a high weight to dense regions, 
and therefore, high weight to small scales. 

4. Spatial locality of the 2MASS galaxies 

4.1. Second order correlations 

We now turn to the major task of this paper - to check whether the number density field 
of the 2MASS galaxies is spatially localized. We start by looking at kj as given by eq.(14). 
The result of K,j vs. 9 is shown in Fig. 3. Here we consider only modes with j = (j, j) and 
j = Q - 10. The error bars are estimated by 90% confidence level for all 28 data squares. 

Fig. 3 shows that the WFC correlation function (ej^ej^i') is substantially different from 
The WFC correlation function is spatially localized for all cases considered. All 
nonlocal (|1 — l'| > 1) correlations are much less than 1, and practically equal to zero within 
the error bars. In other words, the correlation length in terms of the position index 1 is zero. 

The spatial locality is not a surprise on large scales for which the field basically remains 
Gaussian. However, Fig. 3 shows that spatial locality holds on scales as small as = 
28.2/2^° ~ 0.03 degrees. This result is in contrast with Fig. 2, in which the correlation 
function at |1 — l'| = 1 for all j is always comparable with that of |1 — l'| = 0. 

Therefore, it is appropriate to describe the correlation function by (ej lej/ 1/) oc 5/^,. 

One cannot conclude that the spatial locality is due to the dynamical evolution because 
spatial locality can also be given by the Gaussianization due to the central limit theorem 
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(CLT). According to the CLT, a variable given by a superposition of many independent 
identical random variables tends to become Gaussian. If wc consider that a 2-D field is a 
superposition of many identical slides of a 3-D field, the 2-D field variables will be Gaussian- 
izcd. If so, the spatial locality is a trivial result. However, the 2-D WFC ej^i are not subject 
to the CLT. This is because, as discussed in §2.3 ej^ is not a 3-D superposition. This point 
can be directly shown with the non-Gaussianity of variable ej^i. We plot in Figure 4 the 
one-point distributions of the ensemble of ej^i for j = (j, j) with j = 5 - 10 and k — 1...2^ — 1. 
All the distributions of Fig. 4 are non-Gaussian. These one-point distributions generally are 
long tailed, and therefore, the kurtosis of this field is high. 

We also see from Fig. 4 that the non-Gaussianity of the one-point distributions is scale- 
dependent. On large scales j — 5 and 6, the distribution consists of a Gaussian center and a 
power law tail. On scale j — 7, the Gaussian center disappears, and the whole distribution is 
a power law. On small scales j > 8, the distributions still have a power law tail, but the sharp 
center part appears and grows with j. This means that in most places the perturbations 
are small (center part of the one-point distributions), but there are rare events with very 
large perturbations (long tail). That is, there is no Gaussianization of the variable ej^i. The 
spatial locality doesn't arise from the Gaussianization of ej^i. 

In Figure 5, we plot the one-point distributions of the SFCs ej i. These distributions are 
non-Gaussian too. They have a long-tail towards high density, or high value of ej^i. These 
distributions are similar to the one-point distributions given by count-in-ccll. This result 
shows again that both the WFCs and SFCs are highly non-Gaussian. The spatial locality 
cannot be due to the Gaussianization, but a remain of the spatial locality of the initial 
perturbations. 



4.2. Fourth order correlations 

Before looking at the spatial locality of the WFC correlation at higher order we first 
calculate a 4*^ order correlation of the SFC ej i defined by 

\^j,lAej,l+Al/ 

Figure 6 plots the j-dependence of 5?'^ for j = (j, j) with |A1| = 1, 2 and 3 for both the 
2MASS and random samples. The error bars are again the 90% confidence of the 28 samples. 
As expected the random samples always show 5"?'^ = 1, i.e. there is no correlation. On small 

2 2 

scales j — 10, ' for random sample is httle larger than 1 at A = 1. This is because of shot 
noise, i.e. the mean number per cell on small scales is very small. On the other hand, the 
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2MASS galaxy samples show significant non-local correlation S^'^ > 1 for all scales j > 5. 
That is, the correlation length for 4th order statistics of the SFC variables is non-zero. The 
correlation is no localized. On large scales j = 2 - 4, the nonlocal correlation of the 2MASS 
galaxy samples is negligible, because the field is still in Gaussian on such large scales. 

However, the 4*'* order correlation function of the ej^i has completely different behavior 
from the SFC. Similar to eq.(26), we consider a 4th order correlations of the WFCs by 

^2,2 _ \^j,l^j,l+Al/ 
' "(€l)(€l+Al)- ^ ^ 

where Al = 1 — 1'. Eq.(27) is eq.(16) with p = q = 2. Figure 7 presents C?'^ vs. 6 for 
j = (i) i) with j = 6, 7, 8 and 9 for both the 2MASS and random samples. The error bars 
in Fig. 7 are again given by the 90% confidence levels as Fig. 6. Figure 7 shows C?'^ » 1 at 
1 = 1' = for all scales considered. Nevertheless, for all non-local cases |A1| > 0, the 2MASS 
samples behave the same as the random samples within the error bars. Therefore, the 4*^* 
order correlation of the WFCs is spatially localized. Similar to the second order statistics, 
the spatial locality of the 4*^ order WFCs correlation is not due to the Gaussianization of 
Cj 1. At point |A1| = 0, or 6' = 0, the 2MASS samples always have much higher value of C?'^ 
than that of random samples. 

The difference between the WFCs and SFCs 4*^ order correlations can be seen with 
the hierarchical clustering or linked pair approximation. In this approximation, the 4*^^ 
correlations (e|jej i+ai) can be expressed by second order correlations as 

(eJ,ie|,i+Ai) = (^jv)(ejv+Ai) + 2(ej,iej,i+Ai)^ 

+<54[2(e?i)(e?i+Ai)(ej,iej,i+Ai) + (eJ,i)(ej,iej,i+Ai>^ 
+(€i+Ai>(?j,iej,i+Ai>' + 2(ej,iej,i+Ai)'], (28) 

where Q4 is a constant. The first two terms on the r.h.s. of eq.(28) are from unlinked dia- 
grams and the term of Q4 comes from non-linear evolution. Obviously, if the second order 
correlations of the WFCs are spatially locahzed [eq.(14)], the 4*'* and higher order correla- 
tions will also be approximately spatially localized if the linked pair approximation holds. 
Oppositely, the 4th order correlations will not be localized if the second order correlations 
are not spatially localized. 

4.3. Scale-scale correlations of the WFCs 

In the last two subsections, we showed that the gravitational clustering doesn't cause the 
WFC-WFC correlation if the two phase space modes locate in different spatial(x) positions. 
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In other words, to view the dynamical evolution of the cosmic gravitational clustering in 
{k — x) phase space, the gravitational clustering dynamics caused only a very weak coupling 
in the x-direction. Therefore, one can expect, the dynamical coupling should be strong in 
the A;-dircction, which is coupling between modes with the same x, but different k. This is 
the so-called local scale-scale coupling which we measure using 



/~2 ~2 

^2,2 _ \^j,[l/2]gj+l,., , . 

"-^jj+i - /g2 wg2 \ ■ y^^) 

This measures the correlation between modes on scales j and j + 1. The notation [ ] in 
Eq.(29) denotes the integer part of the quantity since the spatial range of the cell (j, [1/2]) 
is the same as the two cells (j + 1, 1) and {j + 1,1 + 1). Therefore, eq.(29) measures the 
scale-scale correlation at the same physical point x. It is the local scale-scale correlation. If 
perturbations on scales j and j -|- 1 are uncorrelated, we have (^[i/2]^+i,i) — (^[i/2])(^j+i,i)' 
and therefore, C'?j+i = 1- 

Figure 8 presents the result of C'?j+i vs. j for the 2MASS galaxies and random sample, in 
which j = (i, j)- This figure shows that the values of C'?j+i for the 2MASS sample are larger 
than 1 on scales j > 6, and also larger than random samples. Therefore, the local scale-scale 
correlation is significant on scales J > 6, or ^ < 0.4 angular degree, corresponding to 1.5 h~^ 
Mpc at the median redshift of the survey. In other words, the gravitational coupling leads 
to the correlation in /c-direction, and the factor in the initial perturbations eq.(lO) is no 
longer hold when nonlinear evolution takes places. 

We can detect nonlocal scale-scale by using 

^2,2 / . ,^ _ (^j,[l/2]gj-+l,l+Al) 
\^j,[l/2]/\^j,l+Al/ 

When I All = 0, Cff_^^{0) is the same eq. (29). When |A1| > 1, C.;f^i{Al) measures nonlo- 
cal scale-scale correlation, because the phase space cell of (j + 1,/ + 1) does not spatially 
overlap with the cell [j, [1/2]). Therefore, the nonlocal scale-scale correlation is measured by 
Cff, i(A/) with I All > 1. 



Figure 9 plots Cjj^j^(A/) vs. 6 for j = (j, j) and j = 6, 7, 8 and 9. The error bars 
are 90% confidence. For all panels, the data points from left to right correspond, one by 
one, to |A1| = 0, 1, 2... The figure shows that the correlation C?j^-^(A/) for the 2MASS 
samples are larger than the random data on the first two data points, but after the third 
data point, the 2MASS samples do not show significant difference from the random data. 
Therefore, the local scale-scale correlations (|AI| = 0) are strong. Non-local scale-scale 
correlations only exist between cells which are spatially nearest neighbors, i.e. |A1| = 1. 
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The scale-scale correlations with |A1| > 1 generally are negligible. For j = 7, 8 and 9, the 
spatial distance between nearest neighbors are only, respectively, 28.3/2-^ = 0.22, 0.11 and 
0.06 angular degrees. On the other hand, the SFCs have strong scale-scale correlation on, 
at least, |A1| = 3 (Fig. 6). In this sense, the scale-scale correlation is also approximately 
spatially localized. 

5. Discussions and Conclusions 

The 2-D distribution of 2MASS galaxies has been studied using the DWT, in which 
the SFC ej 1 is the count-in-cell in physical space, while the WFC ej i is a count-in-cell in 
phase space. We find that the statistical properties of the SFC and WFC variables are 
very different. The former is non-Gaussian and nonlocally correlated, while the later is non- 
Gaussian, but its correlation basically is spatially localized. That is, the sample has the 
following statistical behavior in the phase space 

• the one-point distribution of ej^i is highly non-Gaussian on angular scales less than 1°, 
corresponding to ~ 3.5 h~^ Mpc at the median redshift of the sample. 

• The 2nd and 4th order mode-mode correlations with ej^i are always spatially localized. 
For highest j, the locality holds on angular scales as small as 0.06°, or ~ 0.2 h"-*^ Mpc 
at the median redshift of the sample. 

• the local scale-scale correlations are significant on scales less than 0.4°, corresponding 
to ~ 1.5 h^^ Mpc at the median redshift of the sample. The scale-scale correlation is 
also approximately spatially localized. 

A direct physical meaning of these results is that the cosmic gravitational instabihty causes 
only strong interaction between the modes in different wavebands and in the same spatial 
area, but is weak for modes in different spatial area. 

Because the nonlinear evolution of cosmic gravitational clustering has the memory of 
its initial spatial correlation in the phase space, the observed spatial locahty of the 2MASS 
galaxies provides solid evidence for models assuming that the initial perturbations are spa- 
tially uncorrelated among phase space modes. This result is consistent with the assumption 
that the initial perturbations are Gaussian. Although the Gaussianity of the initial pertur- 
bations on large scales has been extensively tested with the temperature fluctuations of the 
cosmic microwave background radiation, the test given by the spatial locality is effective to 
comoving scales as small as ~ 0.2 h~^ Mpc. 
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The spatial locality has been studied with samples of Ly-alpha forests (Pando, Feng 
& Fang, 2001). Since Ly-alpha forests refer to weakly non-linear clustering, their spatial 
locality can be used to rule out the models of initially non-Gaussian fields which are non- 
spatially localized, but cannot do so for models of initial non-Gaussian which arc spatially 
localized. In this paper, the spatial locality is found for sample referring to fully nonlinear 
regime. The memory of the spatial locality is based on the halo model, for which an initially 
Gaussian field is necessary. In this sense, the spatial locahty of the 2MASS gives stronger 
support to the initially Gaussian assumption than that that of Ly-alpha forests. 

The one-point distribution of 2MASS galaxies (Fig. 3) is very different from that given 
by N-body simulation sample. The later in nonlinear regime generally is lognormal (see Fig. 
6 of Feng & Fang, 2004), while the former is more comphcated (Fig. 3). This difference 
indicates that galaxy distribution is biased with respect to the underlying dark matter field. 
Nevertheless, the galaxy distribution is also highly spatially locahzed. This result indicates 
that relationship between Pg(x) and p(x) should be localized in phase space. That is, the 
stochasticity and nonlocality of the relation between Pg(x) and p(x) arc limited in the cell of 
phase space. For a structure consisting of perturbations on scales k to k -|- Ak, the nonlocal 
size in the relation between Pg{x^) and p(x) should not be larger than |Ax| = 27r/|Ak|. 
Otherwise, the galaxy bias will violate the spatial locality of the galaxy distribution even 
when the underlying field is spatial localized. Hence, from the 2MASS galaxies we can 
conclude that the stochasticity and nonlocality of the bias relation between Pg(x) and p(x) 
probably are no more than that given by the uncertainty relation |Ax||Ak| = 27r. 

The authors thank Dr. J. Pando for his assistance. This research is partially sup- 
ported by the National Basic Research Priorities Programmer of The Ministry of Science 
and Technology of China. 
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Fig. 1. — The number counts dN/dm and completeness C{m) for 1.) the (standard) sample 
of |6| > 30° and 12 < Ks < 13.7 (left panels), and 2.) sample of \b\ > 10° region when 
9 < Kg <14: (right panels). 



0.1 1 

^(degree) 




^(degree) 

Fig. 2. — The SFC correlation functions (pj,iPj,i') vs. 9. The angular distance 9 = 
|1 — r|V800/2-' angular degree, j — 8, 9, 10 and 11, which is the scales •\/800/2 of smoothing 
the data by the scaling function. 
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Fig. 3. — The WFC correlation function Kj vs. e. The angular distance ^ = |1 - l'|y800/2^' 
angular degree, j = (j, j) and j = 5, 6, 7, 8, 9 and 10. 
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Fig. 4. — The one-point distribution of WFCs ej^i for j = and j — 5, 6, 7, 8, 9 and 10. 
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Fig. 5. — The one-point distribution of SFCs ej^i for j = {j,j) and j — 5, 6, 7, 8, 9 and 10. 
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Fig. 6. — The 4*^^ order statistics S-' vs. j with |A/| = 1, 2 and 3. The angular distance 
^ = |1 - I'l ^800/2^' angular de gree. 
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Fig. 7.— The 4*^* order statistics C?'^ vs. 6. The angular distance 6 = |1 - l'|y800/2^' 
angular degree, j = (j, j) and j = 5, 6, 7, 8 and 9. 
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Fig. 8. — The local scale-scale correlations C'?-+i vs. j. The angular scale of j is ^800/2^ 
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Fig. 9. — The nonlocal scale-scale correlations C?j^j^(A/) vs. 9. The angular distance 
e = |Al|v»/2^ j = (j, j) and j = 6, 7, 8, 9. 



